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Abstract 



I discuss the scattering of phonons by a vortex moving with respect to a su- 
perfiuid condensate. This allows us to test the compatibility of the scattering- 
theory derivation of the Iordanskii force with the galilean invariance of the 
underlying fluid dynamics. In order to obtain the correct result we must retain 
0(Vg) terms in the sound-wave equation, and this reinforces the interpreta- 
tion, due to Volovik, of the Iordanskii force as an analogue of the gravitational 

Bohm-Aharonov effect. 
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I. INTRODUCTION 



The problem of computing the transverse force acting on a vortex in a superffuid has 
recently engendered a certain amount of controversy. If the vortex moves at a velocity 
while the superffuid and normal components have asymptotic velocities v s and v n respec- 
tively, then the most general form of the transverse force per unit length that is consistent 
with galilean invariance (i.e. depends only on velocity differences) can be written H 

F = Akz x (v L - v s ) +Bkz x (v L - v„). (1.1) 

Here k is the magnitude of the quantum of circulation about the vortex, 

« = -, (1-2) 
m 

with m the mass of a helium atom, and z a unit tangent to the vortex line in the direction of 
the circulation. In the absence of any normal component, elementary fluid mechanics shows 
that the momentum flux into the vortex core is 

F M = ptotKZ x (v L - v s ), (1.3) 

where p to t is the total mass density. This is the Magnus force . Once a normal component 
is present, however, a variety of different expressions for the coefficients A and B have been 
given in the literature. 

It is generally accepted that the coefficient A is p s [§,|J. An appealing thermodynamic 
argument for this has recently been given by Wexler Q. The controversy stems from the 
claim of Wexler flT], and Thouless et al. [||] that the coefficient of in F, is also equal to p s . 
Since this coefficient is A + B, their claim, if true, would force B to be zero — thus ruling 
out the existence of the second term, the Iordanskii force, which is supposed to originate 
in a left/right asymmetry in the scattering of quasi-particles by the vortex line ||. Sonin 
on the other hand, has presented a detailed review of the scattering of phonons by a 
vortex line at rest with respect to the superflow. His analysis (which has been challenged 
by Wexler and Thouless H[lO|, but which I believe to be correct) shows that the asymmetry 
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arises from a fluid dynamical analog [|TT|] of the Bohm-Aharonov effect ||12|| , and gives the 



coefficient of — v n as p n . Thus B = p n . Combining this value of B with the accepted value 
for A gives the transverse part of the force per unit length as 

F = p s nz x (v L - v s ) + p n KZ x (v L - v„). (1.4) 

This is the most commonly accepted expression for the force. When it is written in this 
form the first term is usually refered to as the superfluid Magnus force and the second term 
as the Iordanskii force. 

Since the total density is p to t = Ps + Pn, equation (|1.4j) may equally well be written 

F = PtotKZ x (v L - V s ) + p n KZ x (v s - v n ). (1.5) 

The first part of ( |1.5| ) is the momentum transfer to the vortex due to the condensate motion 
(and possibly this should be called the superfluid Magnus force), so the second term must 
be the force on the vortex due to phonon scattering. Part of the phonon force is responsible 
for reducing the coefficient of v s from p tot to p s . Notice that this expression for the phonon 
force does not depend on the motion of the vortex line relative to either component of the 
fluid. Although Sonin, who works in the frame = 0, writes p n (y s — v n ) in his expression 
for the phonon force, his analysis of the scattering process is restricted to the situation 
where v s = 0, i.e to the case where there is no relative motion between the vortex line 
and the condensate. The v s part of the force is inferred from the thermodynamic and 
galilean invariance argument given above. A direct demonstration that the phonon force 
is independent of the relative velocity of the vortex and the condensate, and hence that 
coefficient of v s in the phonon force is indeed equal to p n , would provide a useful consistency 



test of the conventional scattering-theory derivation of ( |1.4|) , because the galilean invariance 
that went into deducing this coefficient is not manifest in the linearized sound wave equation. 
In this paper I will provide such a demonstration. In obtaining the result we will find it 
useful to consider the analogy, first pointed out by Volovik ||13|| , of phonon vortex scattering 
with the gravitational Bohm-Aharonov effect where particles are scattered by a spinning 
cosmic string. 



In the next section I will review the acoustic Bohm-Aharonov effect and rederive Sonin's 
results for the phonon force in the case that the vortex is at rest with respect to the con- 
densate. In section three I will extend these results to the case in which the vortex moves 
relative to the condensate, and relate the momentum given to the phonons by the vortex 
to the time delay of signals passing on either side of a cosmic string. In the last section I 
will discuss the apparent conflict between our results and the claims of Thouless et ai, and 
some possible resolutions. 



II. THE ACOUSTIC BOHM-AHARONOV EFFECT 



The scattering of phonons by a superfluid vortex was first studied by Fetter |l3 |. The 
wave equation used by most recent authors fTHKl is 



^VV-2(v.V)f. (2,) 
Here <ft is the velocity potential, v = v„ is the velocity field of the vortex 

((v v ) x , (v v ) y ) = (-2=7-2, TT - 2 ) > ( 2 - 2 ) 

and c the speed of sound. 

When the sound field has harmonic time dependence, <f)(r,t) = e~ lu}t <j)(r), equation 
Ah becomes 



-u 2 (f) = c 2 V 2 (f) + 2iuj(v v - V)0. (2.3) 

We will be interested in effects only to first order in the circulation k, therefore it is natural 
to add harmless 0{y 2 ) terms to ( |2.3| ) so that it becomes the Schrodinger equation for unit 
charge particles minimally coupled to a gauge field A = %v v , viz. 

,1 _ 2 1 — - u 



- oo 2 <f) = c 2 (V + i- 2 v v ) 0. (2.4) 
Notice that this rewriting requires V • = 0. 
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Equation (|2.4| ) describes the Bohm-Aharonov interaction of particles with a threadlike 
tube of magnetic flux in the gauge where V ■ A = 0. The total flux in the tube is 



Here the integration contour surrounds the vortex, which we have taken to be at the origin 
of our coordinate system. We will use the symbol a to denote the ratio of this flux to the 
Dirac flux quantum, $o — 27r. 

In their original paper [|l2j Aharonov and Bohm provided a partial-wave series expansion 
for scattering of a plane wave by the flux tube. Figures |l] and ^| are numerical plots of the 
sum of the first forty terms in this expansion for the cases where a plane wave is incident 
from the right on flux tubes with a equal to 0.25 and 0.5 respectively. These plots should 
be compared to the ripple tank photographs of surface waves interacting with a "bathtub" 
vortex in pT| . The most noticable feature in both the photographs and the plots is the 
"seam" or "tear" in the wavefunction downstream of the flux tube. The incident plane wave 
is cut in two by the flux tube and, apart from diffraction effects, the upper and lower halves 
of the incident wave propagate parallel to each other but with a relative phase shift of a 
fraction a of a wavelength. It is in the region of the seam that the transverse momentum 
imparted to the beam by the flux resides. Indeed in Fig [I] one can plainly see that the 
wavefronts are directed slightly upwards in this region. 

The time- average momentum density in the sound wave is (pmVm) where vm = V(/> is 
the fluid velocity due to the sound wave and 



The angular brackets denote a time average. (See the appendix for a derivation the expres- 
sion for P(i).) For a plane wave 




(2.5) 



|40 + y V0}. 



(2.6) 




(2.7) 



and with the background flow v vanishing, we have 



5 



(P(i)V(i)) = ip (O )^|0o| 2 k. (2.8) 



More generally we find 



(P(i)V(i)) = (0*V0 - (V0* )0) , (2.9) 



where u r is the frequency of the wave relative to the fluid. (Notice that ( |2.9| ) is not the 
"gauge invariant" form of the current for our minimally coupled Schrodinger equation.) 
Once we are out of the region where v is significant, we can write the momentum density 

as 

(P(D V (D> = ^P(o)^l0o| 2 V X , (2-10) 

where <p(x) = \(f)o\e lx ^ r \ If we temporarily ignore the reduction in the amplitude of the sound 
wave in the seam region, we can find the total transverse momentum by integrating the y 
component of this momentum from one side of the seam to the other along a line parallel 
to the y axis. The total transverse momentum per unit length at abscissa x is therefore 

(Pv) = ^(o)^|0o| 2 A X (x), (2.11) 

where Ax{x), the phase difference across the seam, is zero long before the sound wave 
interacts with the vortex, and 



Ax(ar) = 2na = %k (2.12) 
<r 

well after the sound waves have passed the vortex. In this manner the transverse momentum 
is found by examining the wave at large (but not infinite) impact parameter, and the result 
is insensitive to details such as diffraction effects. 

The transverse momentum per unit length of the seam can also be written 

Qp(o)^o| 2 *:) • I ■ ^« = O'ph) • I ■ ^k, (2.13) 

where (j p h) is the mass current, or momentum density in the unperturbed sound wave. If a 
finite pulse of sound is sent past the vortex then a length of seam equal to the group velocity 
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of the waves (here c) is created every second. The transverse force is the rate of creation of 
transverse momentum and this is 

(£±> = 0'ph>'^-^«-c=0'ph>« (2.14) 

in agreement with ref ||. 

The mass flux due to phonons in the two fluid model is 

(j P h) = Pn(v n - v s ), (2.15) 



and we can find the thermal average of the phonon force by substituting this in ( |2.14| ). Since 



we have so far assumed that v s is zero we have, however, only established the p ra ftv n part 
of the Iordanskii force. 

A more rigorous approach to computing the momentum given to the phonons exploits 

the momentum flux tensor II? hon . The relevant terms are 

l j 

n f, h ° n = P(o)( v (i)i v (i)j) + ( v v)i(P{i)V(i)j) + («„)iW) v (i)i>- (2.16) 
The only terms contributing to Tl X y 0n = n^ 011 to first order in v„ turn out to be 

nS on = « on ) Vv=0 + 1 Ux + *(«.),) (UM, (2.17) 

while T[Py 0n is of at least second order in v„ and can be neglected. 

For most of the x, y plane we may use the eikonal approximation for the phase \i 
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X(x,y) = kx-— (v v ) x dx'. (2.18) 



C J — oo 



Here the integral is taken along the line from the point (—oo,y) to (x,y). (The eikonal ap- 
proximation will be described further in the next section.) From the formulae for (v v ) x , {v v ) y 
in Q we find 

/ (v v ) x dx = -— f, tan- 1 r7 + - . 2.19 
J-oo 2ir \y\ \ \y\ 2 J 

This expression is continuous across the negative x axis, but jumps discontinuously by k 
across the positive x axis. We therefore find that the eikonal phase has the expected ujk/c 2 



discontinuity across the seam. The physical wave fronts, of course, smoothly interpolate 
the phase across the seam. We immediately see that outside this interpolating zone, in the 
region where the eikonal approximation to the phase is valid, we have 

dyX+hv v ) y = 0. (2.20) 

This means that U xy is zero outside the seam region. Thus the flux of p y through any curve 
vanishes unless it intersects the seam. Indeed we find that the only regions that have a net 
transverse momentum flux out of them are those that include the vortex. For these the 
momentum flux is entirely due to the interpolating phase and comes out to be ((j p h)x)^ as 
found by the previous, more intuitive, method. 

The transverse momentum flux tensor vanishes because the transverse component of k 
acquired by interaction with the flow is cancelled by the transverse component of v„ when 
they are combined to form the group velocity 

£ = £k + v.. (2.21) 

The classical phonon trajectories therefore remain straight, and, just as in the ordinary 
Bohm-Aharonov effect, the transverse momentum is consequence of the wave-particle dual- 
ity. 

A more detailed analysis would take into account the reduction of the amplitude of the 
sound wave in the seam region. It is well understood from the theory of Bohm-Aharovov 



scattering |15j that the effect of this is to replace $ by sin $ in the force equation. (Inspection 
of Fig ^ where $ = ir shows that that the sound wave amplitude is exactly zero in the seam). 
Using c = 230 ms" 1 for the speed of sound in liquid helium, we find 

a = ^- = 0.035£ phon (2.22) 
c z m 

where -Ephon is the energy of the phonon measured in degrees Kelvin. We see that $ is 
small at temperatures below 0.2 K where phonons dominate the scattering process, so the 
correction will be unimportant. 
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III. MOVING THE VORTEX 



So far we have merely reproduced the results of ||. We now extend our analysis to 
the case in which the vortex moves with respect to the condensate. So as to retain a time 
independent equation we will keep the vortex fixed at the origin, but allow a non-zero 
asymptotic v s . For simplicity of description we will initially consider only the case where 
the uniform v s is in the direction of propagation of the sound wave, which as before we take 
to be the +x direction. We will write v s = C/x. In this case the length of seam created per 
second is c + U. We need to find the phase shift between two halves of the wavefront to 
complete the computation. 

It is tempting to simply replace the in (|2.4j ) by v = C/x + v„, but this will not serve 
to give the correct answer. Because of the Doppler shift, the frequency is now related to the 
wavelength by uj = (c + U)k, so for a wave with the same k the frequency, and hence the 
effective flux 

$ = ^fv v -dr (3.1) 

is increased, and this makes the phase shift larger. This is not what we expect, and is 
incorrect. The terms added to the sound-wave equation to make it into the minimally 
coupled Schrodinger equation are no longer harmless. This is because even when we work 
only to 0{v v ) accuracy, we must not neglect 0{U 2 ) terms. 

We require a more accurate equation. In the appendix we show that the relevent equation 



is that given by Unruh |1|[17] 



(! + H£(! +v ' v V= v " v ^ (3 - 2) 

who interprets his equation as that of a scalar field propagating in a curved space-time 
background. 

We can find the phase offset by working at large impact parameter where <fi is essentially 
independent of y, and where also d x v v is negligeable. We can therefore set 
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= p(x)e ikx - iuJt (3.3) 

with ijj = (c+U)k, and expect <p to be slowly varying. The equation for becomes 

- u 2 <p - 2iuv x d x (f) + v 2 x d 2 xx <p = c 2 dl x (f) (3.4) 

where v x = (v v ) x + U. Since ip(x) is slowly varying we can ignore d xx ip(x). Doing so we find 
that ( |3.4| ) reduces to 



(v v ) x kcp - i(U + c)d x (p = 0. (3.5) 

(I have ignored terms containing v v in the coefficient of d x tp because they will not affect our 
result to O(k).) This gives 

(p{ x ) = e~ ix{x) (3.6) 

with 



\{x) = j^— / (v v ) x dx'. (3.7) 

U + C J-oo 

We see that as x — > +oo the phase offset of the two wavefronts becomes 

The factor of U + c cancels against the length of seam being created per second to give 

(P±) = P h)« (3.9) 

as before. 

This result can be confirmed by examining the momentum flux tensor. The 0(v„) part 
of U xy is now 

= {ic + u)^d yX + (v v ) y ^ (UM, (3.io) 

where we have included a non-zero contribution from U{p(i)V (i) y )- Once again we see that 
outside the seam region the gradient of \ cancels the (v v ) y advective term, and that the dis- 
continuity across the seam provides the momentum flux obtained in the previous paragraph. 
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The wavefront offset can also be calculated from the time delay between phonons passing 
on other side of the vortex. Since the phonon trajectories with large impact parameters are 
hardly deflected, we can find them as the null geodesies of the simplified form of the Unruh 
metric ( [A17p 



ds 2 = — j— (dx — (v + c)dt) (dx — (v — c)dt) — dy 2 — dz 2 ^ . (3-11) 
The null geodesies are given by 



< ^ = v±c. (3.12) 
dt y ' 



In our present case, the time of arrival of a signal at the point (x, y) is 



t(x,y) = I — , r = const. - / (v v ) x dx + 0(v 2 v ). (3.13) 

-oo U+C+(V V ) X (u+cyj-oo 



We convert the time delay to a phase shift by multiplying by u — (U + c)/c. Again we find 
that the relative phase shift between waves that pass above and below the vortex to be 

U + c 



As described by Volovik [IT3J this separation in the time of arrival for signals passing abitrarily 



far from the vortex is characteristic of a "spinning cosmic string" . 

Matters become slightly more complicated when the uniform background superflow v s 
is not oriented parallel (or anti-parallel) to the incident phonon flux. After a little work we 
find that the eikonal equation becomes 

(U 9 - V)x + k- Vt , = 0, (3.15) 

where XJ g = ck + v s . The phase discontinuity seam no longer lies exactly along the x axis, 
but instead points in the direction of U 9 . The rate of transverse momentum production does 
depend on the angle this vector makes with the x axis, but the effects are 0(U 2 /c 2 ). They 
do not seem to be worth working out in detail since at this order we should also include 
compressibility effects and the dependence of p n and p s on (v n — v s ) 2 /c 2 . 
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IV. DISCUSSION 



We have computed the force on a vortex due to the scattering of phonons in the case 
where the vortex is moving with respect to the condensate. We have shown that to first order 
in |v s — v_l|/c the rate of transverse momentum production is independent of the relative 
velocity of the vortex and the condensate, and that the Bohm-Aharonov phase shift of the 
phonon wavefront passing on either side of the vortex leads to a force 

Fphon = PrAZ X (V s - V n ), (4.1) 

which is consistent with galilean invariance. This supports the view that the Bohm-Aharonov 
analysis of the phonon force is correct. 

The expression we have obtained for the phonon force leads to the commonly accepted 
form of the Iordanskii force and so disagrees with the claims of Thouless et al. that it 
vanishes identically. Their claim is based on earlier work by Thouless, Ao and Niu (TAN) 
18] which establishes a general theorem relating the force on a vortex to the circulation of 



momentum at infinity. Since there are no impurities present, this theorem should hold here 
and we have a puzzle that needs to be resolved. 



In response to criticism of their claims p9|j20[| , Thouless et al. |2T[ have given two possible 



explanations for the discrepancy. The first suggests that the computation of the cross section 
asymmetry is mathematically flawed because a conditionally convergent series is mistreated. 
If this were correct, their objection would also apply to the case where = v n , but we have 
shown here that the scattering asymmetry correctly predicts the reduction of the coefficient 
of v s from p tot to p s . This explanation seems unlikely therefore. 

The second possible explanation is more subtle. In the calculations presented here, and in 
all earlier work on phonon scattering, the incident flux of phonons is calculated by assuming 
a thermal phonon distribution derived from the asymptotic v n and v s , i.e. a distribution 
that does not seem to take into account the effect of the local vortex flow field v^. Because 
the assumed phonon flux to the left and right of the vortex is the same, it appears that the 
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phonons make an equal reduction in the fluid momentum on either side of the vortex. In 
other words they appear not to change the value of the total momentum circulation 



so that it remains p tot h/m instead of being reduced to p s h/m. If this were correct there 
would be circulation in the normal fluid as well as in the superfluid component. If we include 
this unphysical circulation in the TAN formula, then it agrees with the calculated Iordanskii 
force. 

The problem with this explanantion is that the motion of the phonons can be derived from 
a hamiltonian, H = c|A;|+v v -k. Consequently Liouville's theorem applies to the distribution 
function. Even in the absence of phonon-phonon interactions, a phonon distribution function 
that is in thermal equilibrium will remain in equilibrium and so phonons arriving from 
infinity will have their local distribution modified by the flow field so as to correctly bring 
about the expected reduction in the circulation. Although there are an equal number number 
of phonons passing to the right or to the left of the vortex, those moving against the superflow 
are going slower, so their number density is higher. This means that the phonon momentum 
opposing the ptot^v is left-right asymmetric as it should be. 

There remains a contradiction, therefore. 

What are we to conclude? The TAN arguments implicitly make use of the change in 
the total momentum of the medium outside the vortex. In fluid mechanics it is well known 
that the total momentum associated with a system of vortices is ill-defined, being given 
by a conditionally convergent integral. This problem is usually dealt with defining the 
impulse [0 of the vortex system. The impulse is defined in terms the velocity potential 
in the vicinity of the vortex system, and does not have contributions from the effects of 
distant boundaries (the ultimate origin of the ill-defined momentum). Perhaps this is at the 
root of the problem. On the other hand the present calculation might be also described as 
naive. We are not distinguishing between the true newtonian momentum and the pseudo- 
momentum |22j possessed by the phonons. However pseudo-momentum is usually exactly 




(4.2) 
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what is needed for computing forces on immersed objects. Clearly, more work is needed to 
resolve the paradox. 
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APPENDIX A: THE UNRUH EQUATION 

The homentropic potential flow of a fluid is derivable from the action 

S = J d*x |p0 + ip(V0) 2 + V(p)} . (Al) 

Here p is the mass density and V0 = v, the fluid velocity. Varying S with respect to gives 
the continuity equation 

p + V-(pV0)=O, (A2) 
while varying with respect to p gives Bernoulli's equation 

+ ^(V0) 2 + Mp) = O, (A3) 

with p(p) = dV/dp. 

In order to consider the propagation of sound waves in the background flow, set 

= 0(0) + 0(1) 

P = P(o) + P(i) + • • • (A4) 

where 0(o) and p(o) obey the equations of motion, and 0(i) and p^ are small amplitude 
perturbations. Expanding S to quadratic order in the perturbations gives 
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S = S + jd 4 x |p ( i)0(i) + X - (j-^j p 2 {1) + ^P ( o)(V0 (1) ) 2 + P(1) v • V0 (1) | . (A5) 

(The terms linear in the perturbations vanish because of our assumption that the zeroth 
order variables obey the equation of motion.) Here v = V( ) = V</>(o). The speed of sound, 
c, is defined by 

c 2 dp 
P(o) dp 

or more familiarly 



(A6) 

P(0) 



^ - %■ < A7 > 

Since the new action is quadratic in p^, we can eliminate it via its equation of motion 

P(i) = -^{0(D + v-V0 (1) }. (A8) 
We find the effective action for the sound waves in the background flow v to be 

S(2) = J d*x {ip ( o)(V0 (1) ) 2 - + v ■ V0 (1) ) 2 } . (A9) 

After changing an overall sign for convenience, we can write this as 

S = J d A x ^ V=0<r^(i)^(i), (A10) 

where 

v^r^ji 1 ' T vT ). (aid 

c \ V, vv J — c z l J 

(We use the convention that greek letters run over all four space-time indices with = t, 
while roman indices refer to the spatial components.) 
The resultant equation of motion 

1 



-9 



duV^ggTd^M = o, (A12) 



is 
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(I + v ■ v ) p -f (I + v ■ v ) = vc ° (0) V0(1)) - (A13) 

This equation and its interpretation as the wave equation for a scalar field propagating in 
the background space-time metric (|A11|) is due to Unruh P^JT! 



In four dimensions we have y/—g = P 2 q)/ c anc ^ 



c 2 — t> 2 , v T 



ft. = ®L I | ■ (AM) 



c 



The associated space-time interval is therefore 

ds 2 = ^ {c 2 dt 2 - 5ij(da? - v i dt)(dx j - v j dt)) . (A15) 



Up to the overall conformal factor ^ we see that c and — v l play the role of the lapse func- 
tion and shift vector appearing in the Arnowitt-Deser-Misner (ADM) formalism of general 
relativity ||23|| . A conformal factor does not affect null geodesies, and so variations in p(o) do 
not influence the ray tracing for the sound waves. 
It is also sometimes convenient to write 

= £ffl { (c > _ ^ (i t ( Sv + J^) ^} . (A16) 

When v is in the x direction only, we can also rewrite ds 2 as 

ds 2 = BR {_ _ ( v + c )dt) (<fc - (v - c)cft) - cfc/ 2 - dz 2 } . (A17) 



This shows that the x — t plane null geodesies coincide with the expected characteristics of 
the wave equation in the background flow. 
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FIGURES 




FIG. 1. Bohm-Aharonov scattering for a = 0.25 
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FIG. 2. Bohm-Aharonov scattering for a = 0.5 



20 



